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ABSTRACT 

We calculate the excitation and dissipation of low-frequency tidal oscillations in uniformly rotating 
solar-type stars. For tidal frequencies smaller than twice the spin frequency, inertial waves are excited 
in the convective envelope and are dissipated by turbulent viscosity. Enhanced dissipation occurs over 
the entire frequency range rather than in a series of very narrow resonant peaks, and is relatively 
insensitive to the effective viscosity. Hough waves are excited at the base of the convective zone 
and propagate into the radiative interior. We calculate the associated dissipation rate under the 
assumption that they do not reflect coherently from the center of the star. Tidal dissipation in a 
model based on the present Sun is significantly enhanced through the inclusion of the Coriolis force 
but may still fall short of that required to explain the circularization of close binary stars. However, 
the dependence of the results on the spin frequency, tidal frequency, and stellar model indicate that 
a more detailed evolutionary study including inertial and Hough waves is required. We also discuss 
the case of higher tidal frequencies appropriate to stars with very close planetary companions. The 
survival of even the closest hot Jupiters can be plausibly explained provided that the Hough waves 
they generate are not damped at the center of the star. We argue that this is the case because the 
tide excited by a hot Jupiter in the present Sun would marginally fail to achieve nonlinearity. As 
conditions at the center of the star evolve, nonlinearity may set in at a critical age, resulting in a 
relatively rapid inspiral of the hot Jupiter. 

Subject headings: binaries: close — planetary systems — stars: oscillations — hydrodynamics — waves 



1. INTRODUCTION 
1.1. Tidal interactions involving solar-type stars 

The tidal interaction between a star and an orbiting 
companion can lead to a significant evolution of the sys- 
tem over astronomical timescales if the orbit is suffi- 
ciently close. When the star experiences a periodically 
varying gravitational potential, an oscillatory tidal dis- 
turbance is generated in the fluid. Dissipation of the 
tide, by whatever mechanism, is directly associated with 
a secular transfer of angular momentum between the spin 
and the orbit as well as a loss of energy from the system. 
All of the spin and orbital parameters therefore evolve 
as a result of tidal dissipation. 

In binary stars the tendency of tidal evolution is usu- 
ally towards a double synchronous state in which the or- 
bit is circ ular and bo th stars spin at the same rate as the 
orbit (e.g. lHutlll981 l). In this state each star experiences 
only a steady tidal distortion and no further dissipation 
or tidal evolution occurs. Observations of stellar popula- 
tions of different ages provide clear evidence for the ongo- 
ing tidal circularization of close binaries, because in the 
older samples circular orbits are found for a wider range 
of orbital periods (Meibom fc Mathieu 2005) . Measure- 
ments of spin are more lim ited but do provide evidence 
of tid al synchronization (Meibom. Mathieu. fc StassunI 
l2006h . 

Many of the known extrasolar planets orbit sufficiently 
close to their host stars to allow significant tidal evo- 
lution. The closest planets, commonly known as "hot 
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Jupiters" or "hot Neptunes", have circular orbits, prob- 
ably mainly as a result of the tidal dissipation within the 
planet rather than that within the star. Although these 
planets are expected to be tidally synchronized, the host 
stars, which are typically similar to the Sun in mass and 
age, are not. The reason is that the orbital moment of 
inertia of a hot Jupiter or Neptune is at most compa- 
rable to the spin moment of inertia of a solar-type star. 
The tidal torque therefore leads to orbital migration ac- 
companied by a modest change in the spin rate of the 
star, which is in any case under the control of magnetic 
braking. 

The direction of orbital migration is away from the 
corotation radius of the star. When the system is young 
and the star rapidly rotating, the inward migration of a 
hot Jupiter driven by its interaction with the protoplan- 
etary dis k may therefore be halted by the tide it r aises in 
the star (|Lin. Bodenheimer. fc Richardson|[l996f ). Once 
the star has spun down, however, tidally driven inspiral 
thre atens the surviv al of hot Jupiters such as OGLE-TR- 
56b (|Sasselovll200l . 

1.2. Timescales of tidal evolution 

The efficiency of tidal dissipation is often 
para metrized by a dimensi onless quality factor Q 
(e.g. iGoldreich fc Soterl Il966f ). which reflects the fact 
that the star undergoes a forced oscillation and dissi- 
pates a small fraction of the associated energy during 
each oscillation period. Since Q always appears in the 
theory in the combination 
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where k2 is the potential Love number of order 2 (a 
dimensionless measure of the interior density profile of 
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the star), it is more convenient to use this combination. 
(For a homogeneous body without rigidity, k2 — 3/2 and 
Q' = Q.) 

The tidal potential experienced by the star can be writ- 
ten as a sum of rigidly rotating components in the form 
of solid spherical harmonics, 

Re['ir%"\e,(j))e-"^'] , (2) 

where (r, 9, (p) are spherical polar coordinates in an in- 
ertial frame of reference centered on the star, uj is the 
frequency in that frame, and 4* is an amplitude. The 
tidal frequency experienced by the star is 

Lo = io — mfi, (3) 
where ft is its spin frequency.^ Q' is a function of £, m, 
and Lo, but can be assumed to be independent of ^! if 
linear theory applies. In practice ^ = 2 is dominant, and 
the azimuthal wavenumber can be restricted to to = 0, 
1, or 2. 

Although the potential applications of tides are very 
broad, we are concerned here with two problems con- 
nected with observations: the circularization of close bi- 
nary stars and the inward migration of hot Jupiters. 

In a binary star with mean motion n > 0, tidal dissi- 
pation leads to an evolution of the semimajor axis a and 
eccentricity e at the rates 
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where we neglect fractional corrections of second order 
in eccentricity and obliquity.'' Here M and R denote 
stellar mass and radius, and Q'i^Q is the value of Q' 
of star i for £ = 2, azimuthal wavenumber m, and tidal 
frequency 6j. For two identical stars that are already 
synchronized (f2i — Q2 = n), eccentricity is damped on 
the circularization timescale 

e 2 / a \ 5 Q' 
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^ We restrict our attention to uniformly rotating stars. By 
quency" we always mean angular frequency. 

^ Under the assumption that 1/Q' oc equation JSj repro- 
duces the result of iDarwirn (II88O,) that eccentricity is excited for 
f2/n > 18/11 and damped otherwise. Note that the discontinuities 
in equations (|4]l and Q do not occur in practice because 1/Q' must 
vanish as <i tends to zero in the case of a uniformly rotating body. 
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Fig. 1. — Circularizatio n period versus age for the s tellar clusters 
and populations listed in IMeibom fc Mathieul II2005I) . The curves 
show the period for which the circularization timescale Q is equal 
to the age for various values of Q' . 

where P is the orbital period and p is the mean density 
(p0 « 1.41 gcm""^). The timescale for synchronization is 
indeed much shorter. Equation assumes a common 
value of Q' for the various tidal components, but this 
assumption may not be so important because the w = 
3n— 2r2 component dominates unless it has a much larger 
value of Q' . 

The relative contribution of the tide in star 2 
is (Qi/Q2)(Mi/M2)2(i?2/i?i)^. For late-type stars 
R cx Af^/^ approximately, giving a ratio of 
{Q'ilQ'2){^h/Mif. The tide in the secondary star is 
therefore less important if its Q' is the same. 

Fig. [T| reproduces the observational data from the 
study of Meibom & Mathieu (2005). Ahhough the trend 
is not particularly clear it is roughly consistent with the 
frequently quoted value Q' w lO*'. 

To determine the rate of inward migration of a hot 
Jupiter, we consider a synchronized planet of mass Mp 
in a circular orbit around a star of mass M^, in which 
case 
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The inspiral time for a planet into a star similar to the 
Sun (assuming Q'^ to be independent of frequency) is 
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The extrasolar planet of shortest orbital period found 
to date is OGLE-TR-56b with P = 1.21 d. It orbits a 
star of mass 1.0 4 Mq whose age is estimated as 3 ± 1 Gyr 
()Sasselovll2003[ ). If the value Q' « 10^ appropriate for 
binary stars is employed in this context, the inspiral time 
is only about 0.036 Gyr. 

The observations raise two basic problems for theorists. 
First, an efhcient dissipation mechanism must be found 
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to explain a Q' as small as 10^ for the binary circulariza- 
tion problem. Second, the same mechanism must not op- 
erate so efficiently in the hot Jupiter problem, unless we 
are extremely fortunate to observe OGLE-TR-56b and 
others like it. 

1.3. Equilibrium and dynamical tides 

The analysis o f tid es in s tars was dev eloped principally 
by IZah^ ([1966303, H^Tfll [TotI UqtI . In the theory of 
the equilibrium tide, the tidal potential induces a time- 
dependent but quasi-hydrostatic bulge. Owing to the 
slow oscillation of the bulge a certain velocity field is re- 
quired within the star. Dissipation occurs in convective 
zones because an effective viscosity associated with tur- 
bulent convection acts on this velocity field. The result- 
ing value of l/Q' is proportional to the effective viscosity, 
which ought to be reduced when the tidal period exceeds 
the convective timescale because then only part of the 
spectrum of turbulent fluctuations damps the tide. 

The assumption of a quasi-hydrostatic bulge appears 
appropriate because the tidal frequency is typi- 
cally much smaller than the dynamical frequency a;* = 
(G'M/i?3)i/2 ^Yie star, which means that fundamental 
and acoustic modes are driven well below their natural 
frequencies. However, stars support other types of wave 
of much lower frequency. The theory of the dynami- 
cal tide determines the excitation of such waves by tidal 
forcing and their dissipation by non-adiabatic processes, 
which can be effective if the wavelengths are short. 

Radiative zones are stably stratified and typically have 
a Brunt-Vaisala frequency N that is comparable to uj^ 
and much greater than Radial motions are strongly 
inhibited when |a)| ^ TV. On the other hand \uj\ is typi- 
cally comparable to 17, so the Coriolis force can be impor- 
tant. In uniformly rotating stars the resulting solutions 
are called Hough waves and have a short radial wave- 
length. These reduce to gravity waves in non-rotating 
stars, as considered in Zahn's work. 

Convective zones are usually adiabatically stratified to 
high accuracy and radial motions are not inhibited. At 
tidal frequencies < 2|f2| in uniformly rotating stars, 
convective zones support pure inertial waves, for which 
the Coriolis force is the restoring force. These have an 
entirely different character from Hough waves, although 
a subset of Hough waves are sometimes called "ine rtial 
modes" (e.s. ISavoniie. Papaloizou. fc Albertslll995f ). In 
a non-rotating star, or for 2\fl\ < \6j\ <C w*, convective 
zones do not support waves. Even in this case, as dis- 
cussed below, the tidal response is not described by the 
equilibrium tide. 

1.4. Previous theoretical approaches 

Several papers have previously presented calculations 
of the efficiency of tidal dissipation in solar-type stars 
under a variety of assumptions. As for the theory of 
the equilibrium tide, some controversy has surrounded 
the correct way to reduce the effective viscosity when 
the tidal per iod exceeds the conve ctive timescale, as 
discussed by iGoodman fc OhI (|1997D . The theoretical 
arguments generally suppor t, at least in an a p proxi - 
mate sense, the formula of Go ldreich fc Keeleyl (|1977l ) 
[cf. equation pH)) below], which drastically reduces the 
effective viscosity at high frequency. Recently, however. 



Penev et al.l (|2007[ ) find some support for the formula of 
ZahnI ( r966bl ) from an analysis of numerical simulations 
of turbulent convecti on. Ev en with the more optimistic 
assumption of Zahi] (|1966b[ ). Q' is too large to explain 
binary circularization. 

Bot h pTerquem et al.l (|1998l ) and IGoodman fc DicksonI 
(|1998f ) extended the theory of the dynamical tide, previ- 
ously developed by Zahn for early-type stars with con- 
vective cores, to solar-type stars with convective en- 
velopes. In this approach, the tidal forcing excites an ir- 
rotational flow in the convective zone an d low-frequency 
gravit y waves in the radiative interior. iTerquem et al.l 
()1998D calculated adiabatic forced gravity waves, which 
propagate through the entire radiative zone, and deter- 
mined the dissipation rate by including non-adiabatic 
effects as a perturbation. Neither radiative damping 
nor turbulent viscosity is very effective in damping the 
waves because they are evanescent in the convective 
zone. The outcome was that tidal dissipation is ef- 
ficient only when the tidal frequency closely matches 
that of a global g mode, but the system would evolve 
very slowly between these resonances. When the dissi- 
pation rate is averaged in time, using an appropriately 
weighted local average in frequency, the resulting value 
is dominated by the behavior off resonance and corre- 
sponds to Q' w 5.1 X 10^(Ptido/10d)-" i5 [based on their 
equation (37)] where Ptide = 2tt/\lj\ is the tidal pe- 
riod. They us ed an effective visco s ity of the form rec- 
ommended by I Goldreich fc KeelevI ()1977[ ) but which is 
in fact larger by a factor of Air^jn the high-frequency 
l imit and closer to that used by Goldreich & Nicholso^ 
()l977t ).^ This value of Q' is substantially larger than 
that calculated on the basis of the equilibrium tide, 
Q' w 1.2 X 10^(Ptido/10d)° °s [based on their equa- 
tion (36)], as that approximation gives an inaccurate de- 
s cription of the disturb ance in the convective zone. 

IGoodman fc DicksonI (jl998l ) did not calculate global 
g modes but considered the possibility that the gravity 
waves would reflect imperfectly from the center of the 
star, where the wavelength is very short and nonlinear- 
ity can occur because the wave energy is concentrated 
into an extremely small volume. Imperfect reflection 
could mean either that the wave reflects from the cen- 
ter with a perturbed phase, or that the wave is substan- 
tially dissipated at that location. In either case global 
g-mode resonances do not occur and tidal dissipation oc- 
curs over a continuous range of tidal frequencies. The 
result is Q' « 1.2 x 109(Ptide/10 d)^/^ [^ased on their 
equation (13)]. (When they converted this result into a 
circularization timescale they apparently underestimated 
it by a factor of about 20.) The stronger frequency de- 
pendence in this expression means that this mechanism 
could be very important for short-period systems. 

Neither of these studies included the effect of the Cori- 
olis force on the tidal disturbance, which might be very 
important if the tidal frequency is not much larger than 
the spin frequency of the star. Incorporating the Coriolis 
force is difficult because it leads to wave equations that 
are not separable in the spherical polar coordinates r 

^ Such factors are certainly debatable as neither of the cited 
papers goes beyond simple dimensional est imates. Neither includes 
the factor of 1/3 appearing in equation mOI I and used by subsequent 
authors. 
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and 6. This separability is very useful because the wave- 
length of the gravity waves can be so short that a two- 
dimensional calculation with adequate resolution is un- 
feasible. (The third, azimuthal, direct ion can always be 
treated by a separation of v ariables.) ISavonije fc Witti 
11002) and IWitte fc Savonije (2002) were able to retain 
separability by using the "traditional approximation" in 
which the latitudinal component —17 sin of the angu- 
lar velocity is neglected in calculating the Coriolis force. 
Indeed, the traditional approximation is justified in the 
radiative region because radial motions are suppressed 
for IcDl <C and t h e rele vant s olutions are Hough waves . 
iSavoniie fc Wifri^ ([200I ) and IWitte fc Savcmiii (pOOl ) 
developed a highly sophisticated model including the ef- 
fects of stellar evolution, magnetic bra king, an d reso - 
nance locking. The res ults of .Savoniie fc Wittd (|2002D 
are similar to those of iTerguem et al.l (|1998D off reso- 
nance, as they adopt the same effective viscosity, but 
inclusion of the Coriolis force allows a richer spectrum of 
resonances with various types of global Hough modes. 

1.5. Purpose of this calculation 

Unfortunately the traditional approximation is not 
valid in the convective zone, which is adiab atically strati- 
fied t o high accuracy. In a previous paper (jOgilvie fc LinI 
12004 hereafter Paper I) we studied the excitation and 
dissipation of tidal disturbances in rotating giant plan- 
ets, which are mostly or fully convective. As noted 
above, for tidal frequencies < 2\Q\, the convective 
zone supports pure inertial waves, which are not cor- 
rectly described by the traditional approximation. The 
importance of inertial waves for tidal dissipation i n giant 
planets has also been recognized by IWul ()2005ai rbl) and 
iPapaloizou fc Ivanovl (|2005). 

Inertial waves have remarkable properties, some of 
which have been elucidated only since the advent of high- 
resolution numerical calculations. They propagate along 
characteristic rays that are inclined to the rotation axis 
at a certain angle, depending on the wave frequency. This 
angle is necessarily preserved in reflections of the waves 
from boundaries, so that a beam is typically either fo- 
cused or defocused in such a reflection. When propagat- 
ing within a spherical annulus (rather than a full sphere) , 
the waves are generically focused on to "wave attractors" 
where intense dissipation occurs fRicutord fc Valdettarol 
[1997; Ricutord, Gcorgcot, & Valdcttaro 2001^)^ It was 
demonstrated mathematically bv lOgilviel (2005) that the 
dissipation rate associated with wave attractors in a 
slightly simplified wave equation is independent of viscos- 
ity in the limit of very small Ekman number. In Paper I 
we found evidence for a similar behavior in the limit of 
small viscosity, although the numerical solutions indicate 
that the dissipation is typically concentrated along the 
rays that emanate from the critical latitude on the inner 
boundary. 

The purpose of the present calculation is to determine 
the excitation of inertial waves in the convective zones of 
solar-type stars and to assess their role in tidal dissipa- 
tion. This occurs both through the damping of the waves 
in situ by turbulent viscosity and also through their ef- 
fect on the excitation of Hough waves in the radiative 
zone. We aim to determine what range of values of Q' 
can be obtained for the binary circularization problem 
and the hot Jupiter problem. 



In some sense our calcula tion is complementary to that 
of lSavoniie fc Wittd (pOOl) . They gave an accurate treat- 
ment of the radiative zone but not the convective zone, 
and assumed that Hough waves reflect coherently from 
the center of the star. We aim to treat the convective 
zone more accurately, and assume that the Hough waves 
do not reflect coherently from the center. 

2. NUMERICAL ANALYSIS 

2.1. Stellar model 

For an accurate and detai led model of th e prese nt Sun, 
we ad opt "model S" from IChristensen-D alsgaar d^ al.l 
In the mixing-length theory of Bohm- Vitensd 



1996). 



19581 ). the convective energy flux density is related to 



the convective velocity v by 
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where 6 — —{d\np/d\nT)p and £ — aHp is the mix- 
ing l ength , = p/pg being the pressure scaleheight. 
[See IZahnI (|1989D for a discussion of the dimensionless 
numerical prefactors.] The convective timescale is then 
r = £/v. 

The turbulent viscosity based on mixing-length theory 
is ^v£, but this should be reduced when the tidal period 
exceeds the convective timescale. We take the effective 
turbulent viscosity to be 
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which is a sm oothed 
iGoldreich fc KeelevI 
form of equation (flOl 
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In fact, 
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functional 



is suggested by a Maxwellian 
viscoelastic model in which r is identified with the 
relaxation time of the turbulent stress. 

Fig. [5] shows the viscosity profile in the Sun for a tidal 
period of 10 d. This period is chosen as being repre- 
sentative of the binary circularization problem. The vis- 
cosity is substantially reduced below the mixing-length 
value in most of the convective zone. For reasonable 
spin periods, the dimensionless Ekman number i'/2ilR^ 
is small but not extremely small (typically of the order 
of 10~^ — lO""*). This means that inertial waves will 
be modestly affected by viscosity, and the problem is 
conveniently accessible to computation with achievable 
numerical resolution. 

Although the interior angular velocity profile of the 
Sun has been measured, we assume the star to be uni- 
formly rotating. This is partly for technical reasons, as 
our numerical method can cope easily only with "shel- 
lular" rotation profiles. In any case it is not clear how 
the rotation profile would change if the Sun were tidally 
synchronized at a different rate in a binary system. 

2.2. Inertial waves in the convective zone 

We solve the forced linearized equations in the con- 
vective zone as in Paper I. The response consists of an 
equilibrium tide plus a dynamical tide. The former sat- 
isfies a second-order ordinary differential equation in r 
which we solve numerically. The latter satisfies partial 
differential equations in r and 9 which can be simplified 
in the low- frequency limit, the simplifications amounting 
to the Cowling and anelastic approximations. The veloc- 
ity field is separated into spheroidal and toroidal parts 
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0.7 0.8 0.9 1.0 

Fig. 2. — Effective viscosity profile in the convective zone of 
the present Sun. The mixing- length value v = ^vt is shown as a 
dashed line, while the value reduced according to equation mOI I for 
a tidal period of 10 d is shown as a solid line. For a spin period of 
10 d a kinematic viscosity oi u = lO^^cm^s"^ corresponds to an 
Ekman number of Ek = u/2Q.B? 1.4 X lO'^. 

and the equations are converted into a large algebraic 
system by Galerkin projection on to normalized associ- 
ated Legendre polynomials (i.e. spherical harmonics) in 
9 and Chebyshev collocation in r. The linear system is 
solved by a standard method for block tridiagonal matri- 
ces, and the total viscous dissipation rate is calculated. 
We adopt stress-free impermeable boundary conditions 
on the dynamical tide at the upper and lower limits of 
the convective zone, as explained in Paper I. 

2.3. Hough waves in the radiative zone 

We determine algebraically the excitation of Hough 
waves in the radiative zone as in Paper I. They are ex- 
cited at the interface between the radiative and convec- 
tive zones, partly by tidal forcing in that region and 
partly by the pressure of the inertial waves acting at 
the interface. FoUowing lGoodman fc DicksonI (1998), the 
model assumes that N"^ vanishes in the convective zone 
and increases initially linearly with distance into the ra- 
diative zone. We assume that the waves are not reflected 
coherently from the center of the star and calculate the 
resulting energy flux. This is converted into a dissipation 
rate as described in Paper I. The behavior of waves near 
the center of the star in the hmit |a)| » is discussed 
in Append ix El where, in a sligh t refin ement of the cal- 
culation of iGoodman fc DicksonI ()1998[ ). we estimate the 
conditions under which the waves become nonlinear. 

3. RESULTS 

3.1. Inertial waves 

In the left half of Fig. [3] we plot, as a function of the 
tidal frequency, the value of Q' resulting from the viscous 
dissipation of inertial waves in the convective zone of the 
Sun with a spin period of 10 d. Only the most impor- 
tant azimuthal wavenumber m = 2 is considered. Also 
shown for comparison is the viscous dissipation rate of 
the irrotational disturbance generated in the convective 
zone when the Coriolis force is omitted. For frequen- 



cies \Cj\ < 2|il| the dissipation is greatly enhanced by the 
excitation of inertial waves, but outside this range the 
Coriolis force has little net effect. 

Where the dissipation rate is significantly enhanced, it 
is relatively insensitive to the viscosity, as found in Pa- 
per I. Elsewhere the dissipation rate is directly propor- 
tional to the viscosity. The value of Q' obtained when 
th e Coriolis force is om itted is larger than that found 
bv lTerquem et al.l (|1998l ). This is partly because of the 
contribution of g-mode resonances to their average value, 
but mainly because their viscosity is larger than ours by 
a factor of approximately 47r^. 

The numerical convergence of these results was verified 
by repeating the calculation with double the resolution 
in each direction. It was found adequate in most cases to 
truncate the Legendre and Chebyshev polynomial bases 
at an order of 100. 

3.2. Hough waves 

Also shown in Fig.[3]is the same information for Hough 
waves excited at the interface between the radiative and 
convective zones. Inclusion of the Coriolis force can ei- 
ther increase or decrease the dissipation rate. 

When the Coriolis force is neglected, our results should 
agree with those of IGoodman fc DicksonI (|1998( ). We 
find their numerical parameter dc to equal —1.19 rather 
than —0.64, which we verified by numerically integrat- 
ing their equation (3). Altogether we find Q' « 3.9 x 
10*(Ftidc/10d)S/3 when spin is neglected. We also ob- 
tain their equation (15) but with {Mi + Al2)/Mi raised 
to the power —5/3 instead of -1-11/6. This discrepancy 
suggests that they may have overestimated the circular- 
ization rate by a factor of 2^/^. 

3.3. Higher tidal frequencies 

In Fig. |4| we plot similar results over a wider range 
of tidal frequencies. As expected, the Coriolis force 
is unimportant whe n \lu\ ^ \n\. As found by 
IGoodman fc DicksonI ([1998), Hough waves (or gravity 
waves, in this regime) provide efficient tidal dissipation 
at high frequencies if they do not reflect coherently from 
the center of the star. 

3.4. Dependence on spin period 

In Fig. [Sjwe present the results for the Sun with a spin 
period of 30 d. Three times the range of uj/il is plot- 
ted, corresponding to the same range of a) as in Fig. |3l 
As expected, outside the range < 2|n|, Q' is almost 
independent of the spin frequency SI at a given tidal fre- 
quency a), because the effect of the Coriolis force is weak. 
Within the frequency range of inertial waves, a smaller 
Q' is generally achieved in the more rapidly rotating star. 
Fig. |6| shows the results for a spin period of 3 d. An even 
smaller Q' results from inertial waves in this case. The 
scaling is approximately 1/Q' oc fl^ for a fixed value of 
uj/n, as suggested by the factor fg defined in Paper I. 
Increasing fl at fixed a)/r2 also decreases the Ekman num- 
ber, resulting in a more structured graph as the system 
is less controlled by viscosity. 

3.5. Lower-mass stars 

We also investigated models of lower mass, as the 
secondaries or even the primaries of "solar-type bina- 
ries" may have masses significantly less than 1 Mq . The 
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Fig. 3. — Dissipation rate, parametrized as a value of Q', as a function of tidal frequency. These results refer to tidal forcing by the 
t = m = 2 solid harmonic in a model of the Sun with a spin period of 10 d. Left: Q' from the viscous dissipation of inertial waves in the 
convective zone. Right: Q' from the excitation of Hough waves in the radiative zone. The dashed lines show the effect of omitting the 
Coriolis force. In this figure only, the dotted lines show the result of artificially increasing the effective viscosity by a factor of 10. 
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Fig. 4. — Expansion of Fig.[3]to a wider range of tidal frequencies. 

greater depth and density of the convective envelopes of 
these stars allows for a more efficient excitation of iner- 
tial waves, in the sense that the dimensionless dissipation 
rate -Dvisc/C^visc (as defined in Paper I) is generally larger. 
However, this effect is offset by the fact that the conver- 
sion factor /q for these stars is smaller as a result of their 
greater mean density. 

In Fig. [7] we plot the values of Q' resulting from inertial 
and Hough waves in a stellar model of mass 0.5 Mq, age 
5 Gyr, and spin period 10 d. Although the details are 
different, the qualitative behavior and the range of values 
of Q' obtained are similar to the case of the Sun. 

3.6. Summary 

In the frequency range \uj\ < 2|f2|, tidal dissipation in 
the convective zones of solar-type stars is substantially 



enhanced through the excitation of inertial waves when 
the full Coriolis force is included. Here Q' can be de- 
creased by up to four orders of magnitude and has a 
complicated dependence on tidal frequency. Typically 
Q' (X ri~^ for a fixed ratio aj/fl. Values as small as 10^ 
can be achieved if the star spins more rapidly than the 
Sun. 

If the Hough waves excited at the interface between the 
radiative and convective zones do not reflect coherently 
from the center of the star as a result of their nonlin- 
earity, they provide another means of dissipation at all 
frequencies. The resulting value of Q' scales generally 
with although this is modified, especially in the 

range \uj\ < 2\^l\, by the inclusion of the Coriolis force. 
Values as small as 10® can be achieved if the tidal period 
is as short as 1 d. The estimates in Appendix |X] sug- 
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Fig. 5. — As for Fig. |3] but for a spin period of 30 d. Note the wider range of li/H, corresponding to the same range of uj as in Fig. |3] 
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Fig. 6. — As for Fig. |3] but for a spin period of 3 d. Note the different vertical scale. 



gest that, in the case of the present Sun, Hough waves 
become nonlinear in eccentric binaries but probably not 
in the hosts of hot Jupiters. Nonlinearity is less likely in 
stars that are younger or less massive than the Sun, but 
is more likely in older stars. 

4. COMPARISON WITH OBSERVED SYSTEMS 
4.1. Tidally induced orbital migration 
4.1.1. Close-in gas giant planets around G dwarfs 

One immediate application of the present analysis is to 
the orbital migration of cl ose-in ext r asola r planets. In his 
study of OGLE-TR-56b, ISasselovl ()2003f ) computed the 
orbital decay timescale a/\a\ (i) by extrapolating a for- 
mula calibrated on the circularization of binary stars, (ii) 
from equilibrium tidal models calculated using the pre- 
scription of Zahn. Q,966b) for the effective viscosity in the 
convective envelope, and (iii) by a similar method using 



the prescription of iGoldreich &: KeelevI ()1977t ). (He set 
the tidal period equal to the orbital period, thereby un- 
derestimating the forcing frequency by a factor of approx- 
imately 2. The tidal frequency associated with the ap- 
propriate m — 2 potential component in the hot Jupiter 
problem is in fact uj — 2n — 2n.) He showed that, while 
the theoretically derived a/\a\ exceeds 10 Gyr with pre- 
scription (ii) and 4000 Gyr with prescription (iii), is it 
less than 1 Gyr with prescription (i) . This discrepancy 
is similar to that found in equilibrium tidal models for 
solar-type binary stars (see Section 14.2. 1|) . 

To the extent that the mechanism and efficiency of 
tidal dissipation remain uncertain, p rescription (i) sh ould 
provide a useful empirical estimate. ISasselovl (|2003l ) also 
estimated the age of OGLE-TR-56 as r* « 3 ± 1 Gyr, 
which is considerably longer than his estimate for a/\a\ 
based on prescription (i). Note further that, owing to 
the accelerating nature of inward orbital migration, the 
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Fig. 7. — As for Fig. |3] but for a star of mass 0.5 Mq and age 5 Gyr. 

inspiral time Tq is considerably shorter than the present 
value of a/\a\ [see equation ©]. While OGLE-TR-56b 
was the only known extrasolar planet with a period close 
to 1 d, a possible interpretation was that it does indeed 
have Ta ^ and we observe it with a small prob- 
ability (~ Ta/Tsf) while it evolves through this period 
range with a rapidly decreasing Ta (oc pi3/3 ^j^jjgp ^j^g 
assumption of constant Q'). However, two additional ex- 
trasolar planets, 0GLE-TR-113b and OGLE-TR-132b, 
have been found with only slightly longer orbital periods 
than OGLE-TR-56b. The values of Ta obtained for these 
planets under prescription (i) are also comparable to or 
shorter than their estimated ages. The period distribu- 
tion of the known planets does not statistically indicate 
that these very short-period planets are at the endpoints 
of their tidal orbital evolution. An alternative interpre- 
tation of these data is that the relevant value of Q' of the 
solar-type host stars may be significantly larger than the 
one that applies to binary circularization. 

Based on the results presented in Section [31 we ex- 
plore the possibility that this dichotomy in the stellar 
Q' may be caused by the differences in the tidal forcing 
frequencies and spin frequencies for the hot Jupiter and 
binary circularization problems. The shortest orbital pe- 
riods of known extrasolar planets are little longer than 
1 d, which is an order of magnitude smaller than the 
circularization periods of mature solar-type binary stars. 
Owing to selection effects, the spin periods of the known 
planet-bearing stars are comparable to or longer than 
that of the Sun. (The absence of significant chromo- 
spheric activity, which is generally correlated with rapid 
stellar spin, has been an important criterion in the se- 
lection of target stars for radial- velocity surveys.) In 
contrast, circularized binary stars appear to attain spin 
synchronization so that their spin p eriod is roughly 10 d 
(Meibom. Mathieu. fc Stassunll2006[ ) . This state is main- 
tained despite the stars' loss of spin angular momentum 
through magnetic braking because this loss can be com- 
pensated by the tidal transfer from the orbit in the case of 
two massive companions. Extrasolar planets have much 
less orbital angular momentum and the tidal interaction 



between the star and the planet cannot compensate for 
the loss of spin angular momentum through magnetic 
braking. 

The important difference is that the tidal forcing fre- 
quency falls within the range of inertial waves in the bi- 
nary circularization problem and well outside it for the 
hot Jupiter problem. In a synchronized binary the eccen- 
tricity tides have frequencies u) = ±17 (cf. Section II. 2p 
and can efficiently excite inertial waves, leading to an 
enhanced tidal dissipation. In contrast, planetary com- 
panions in circular orbits are unable to excite inertial 
waves in their host stars unless $7* > n/2. In Fig. [5] the 
tidal forcing by a planet in a 1 d orbit falls at cD/fi = 18 
and results in Q' 8.9 x 10^" from viscous dissipation in 
the convective zone; this conclusion depends only weakly 
on the stellar spin for \VL\ ^ n. On the other hand. 
Fig. [¥] shows that Hough waves would be very impor- 
tant in the tidal evolution of short-period planets, lead- 
ing to Q' « 1.6 X 10^ for P = Id, provided that they 
fail to reflect coherently from the center of the star. The 
nonlinearity parameter is estimated in Appendix |^ as 
A « 320(Mp/M*)(P/d)^/^ in the case of the present Sun, 
which is about 0.3 for OGLE-TR-56b and somewhat less 
than unity for hot Jupiters generally. The nonlinearity 
increases with stellar mass and age. This estimate shows 
how important it is to determine the fate of Hough waves 
(or gravity waves, in this regime) approaching the center 
for various levels of nonlinearity. 

In contrast, for binary stars in a nearly synchronous 
state in which w w w n for the dominant eccentric- 
ity tide, the Coriolis force strengthens the energy dissi- 
pation associated with both inertial and Hough waves. 
Over most of the frequency range < 2\il,\, the val- 
ues of Q' associated with inertial waves are 2-3 orders 
of magnitude smaller than those outside the range. This 
frequency dependence makes it possible for the dissipa- 
tion of inertial waves to speed up the circularization of 
binary stars (see Section 14.2. ip without significantly in- 
fluencing the orbital evolution of the close-in planets. 

Planets are formed in protostellar disks. Classi- 
cal T Tauri stars are observed to rotate with peri- 



Tidal dissipation in rotating solar-type stars 



9 



ods of a few days (|Bouvier at al.1 [l993l ) and they spin 
up as they contract on to the zero-age main sequence. 
These rapid spin rates are sustain ed for r* < 100 Myr 
(|Bouvier. Forestini. &: AUainl Il997l ). During this early 
phase of stellar evolution, the host stars of the very short- 
period planets may have had fi, > n/2, allowing the ex- 
citation of inertial waves. A relatively small Q' can occur 
in this regime owing to the rapid spin (see Fig. [6]). In 
addition, the properties of the convection may have been 
different. Nevertheless, this phase may contribute little 
to orbital evolution as it is short-lived. 

4.1.2. Close-m brown dwarf around an M dwarf 

A precision radial velocity survey also led to the dis- 
covery of an Mp sin i = 19 Mj planet HD 41004Bb on a 
1.33 d orbit around a 0.4 M© M- dwarf star (jZucker et al.l 
|2003[) . This period is the shortest known for any brown 
dwarf companion, and the estimated age of the system 
is = 1.6 ± 0.8 Gyr. Assuming that the M dwarf fails 
to spin synchronously because of the effects of magnetic 
braking, we find from equation ([8]) that Ta ~ t* when 
Q' w 3.7 X lO"^. Alternatively, if we apply the frequency- 
independent Q'- value (~ 10^), calibrated on the circular- 
ization of solar-type binary stars, the extrapolated value 
of Ta would be nearly two orders of magnitude smaller 
than the estimated age. 

With its low mass, the host star HD 41004B has a 
smaller radiative core than th e Sun. S ince it has no 
measurable rotation speed (But ler et al.l [2006). the tide 
raised by the brown dwarf companion is likely to be well 
outside the frequency range of inertial waves where Q' is 
relatively small. The results in Section [3] indicate that if 
the spin period of HD 41004B is longer than about 3 d, 
viscous dissipation in the convective zone provides a Q' 
much larger than 10^ and poses no threat to the survival 
of the brown dwarf. However, if during the early phases 
of HD 41004B's evolution O was substantially larger than 
its present value, the tidal frequency could have fallen in 
the range of inertial waves. Provided the resulting low 
Q' did not last more than 10* yr, the amount of orbital 
decay would still have been negligible. 

Hough waves can also be excited in the radiative zone 
of HD 41004B. Despite the relatively large mass ratio 
Mp/M^ of this system, the small value of dN/dr at the 
center of the 0.4 Mq star means that the Hough waves 
are unlikely to become nonlinear (see Appendix [A|. This 
is important because if the Hough waves were effectively 
dissipated they would lead to Q' < 10^ and seriously 
threaten the survival of the brown dwarf. 

4.2. Tidal circularization 
4.2.1. Solar-type binary stars 

In the introduction, we indicated that the dissipation 
of tidal disturbances in both the convective and radiative 
zones of non-rotating solar-type stars fails, by about two 
orders of magnitude, to match that needed to account for 
the variation of the circularization period of binary stars 
as a function of their age (see Fig.[T]). These discrepancies 
are particularly noteworthy since the previously obtained 
values of Q' might be underestimated as a result of op- 
tim istic assumptions. For exam ple, in the application of 
the iGoldreich fc KeelevI (Il977f) formula, a factor of iir^ 
was introduced when lTerguem et all (|1998f) replaced d)r 



with r/Ptido- The discrepancies would be still larger if 
the pr escription (fTU|) were applied. Similarly, iSasselovl 
()2003[) set the tidal frequency equal to the orbital fre- 
quency and thereby underestimat ed the suppression of 
the turbulent viscosity. Finally, iGoodman fc DicksonI 
(i l99S) assumed in one section of their paper that Hough 
waves are unable to reflect coherently from the center 
and form global standing modes, thereby allowing effec- 
tive dissipation at all frequencies. 

The results presented in this paper show that the Cori- 
olis force can introduce a frequency dependence in and re- 
duce the magnitude of Q'. The timescale for the stars to 
synchronize is generally much shorter than that required 
for circularization. The tidal frequency Wc = 3n — 2Q 
of the principal tidal component involved in circular- 
ization is larger than that involved in synchronization 
{(lis = 2n — 2n,). Eccentricity can be excited for suffi- 
ciently large spin frequency when lOc becomes nega- 
tive. However, during the spin-down of the star towards 
synchronization, ujs tends to zero, ujc becomes positive 
and tends to Q, and the eccentricity of the binary is 
damped. For this tidal frequency, the results in Fig. [3] 
indicate that the values of Q' associated with inertial 
and Hough waves are 5.7 x 10^ and 1.1 x 10*, respec- 
tively. These values are an order of magnitude larger 
than that obtained for equilibrium and dynamical tides 
in non-rotatin g stars if the gravity w aves reflect from the 
stellar center (jTerguem et al.|[l998[) . But our values of 
Q' are still an order of magnitude too large to account 
for the observational data in Fig. [TJ Thus the incorpora- 
tion of the Coriolis force reduces but does not eliminate 
the discrepancies between theory and observations. 

There are several potential mechani sms that may con- 
tribut e to the circularization process. IWitte k Savonijd 
(|2002^ made an attempt to take the effects of stellar evo- 
lution into account. They suggested that binary systems 
may undergo a complicated evolution through resonant 
tidal interactions. In principle a strong tidal dissipation 
efficiency at a)c ~ could lead to effective eccentric- 
ity damping after the binary stars are synchronized. But 
the protracted maintenance of this feature is needed over 
timescales comparable to Tg. 

The calculations described in this paper are carried 
out under the assumption of uniform rotation. Differ- 
ential rotation is directly observed on the surface on the 
Sun and measured with helioseismological analysis in the 
interior. A preliminary study of models with a modest 
amount of shellular differential rotation indicates that 
the details of the graph of Q'{uj) change significantly 
when differential rotation is introduced, while the typ- 
ical range of values of Q' obtained is similar. 

It is possible that a detailed study taking into account 
the effects of stellar evolution (including spin evolution) 
and differential rotation may be required to account for 
the tidal evolution of binary stars. 

4.2.2. Low-mass binary stars 

Most binary stars do not have equal-mass companions. 
In the solar system, the satellites of gas giants gener- 
ally dominate the circularization process because they 
have much lower values of Q' th an their host planets and 
are more intensely perturbed (|Goldreich fc Soteil 119661 : 
iMurrav fc Dermottlil999) . We now consider the contri- 
bution of a lower-mass companion on the circularization 
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process for binary stars. 

We have computed the rates of tidal dissipation via in- 
ertial and Hough waves for a 0.5 Mq main-sequence star 
of age 5 Gyr (Fig. [7]). This star is less centraUy condensed 
than a 1 Mq star and the fractional radius occupied by 
its radiative zone is smaller. The contribution of Hough 
waves is generally smaller than that of inertial waves un- 
less the tidal forcing frequencies are outside the range 
of inertial waves. Owing to the different internal struc- 
ture, the magnitude and frequency-dependence of Q' in 
this low-mass stellar model differ from that of a 1 Mq 
star. However, when we take into account the mass- 
radius relation, the theoretically determined circulariza- 
tion timescale, in the frequency range of inertial waves, 
is a generally decreasing function of the stellar mass. 

The results of these models can be directly compared 
with some observed systems. The dataset used to derive 
a cir cularization p eriod for the Hyades inclu des Johnson 
311 (jGriffin et al.lfl985l : IMeibom fc Mathieull2005| ). This 
system has a period of 8.5 d and a negligible eccentricity. 
The mass of both of its two components is 0.5 Mq which 
make it an ideal system for a case study of low-mass stars. 
From equation ([6]), we find that each star must have Q' < 
1.8 X 10^ in order for the circularization timescale to be 
less than the age of the Hyades cluster (600 Myr). This 
observationally calibrated value of Q' is more than two 
orders of magnitude smaller than that computed for the 
efficiency of tidal dissipation through both inertial and 
Hough waves. The discrepancy here is larger than that 
found for 1 Mq stars. However, IZahn fc Bouchetl (|1989l ) 
argued that the circularization of these systems can occur 
during the Hayashi phase. 

4.2.3. Close-in brown dwarfs and gas giant planets 

All the theoretical models presented here are for uni- 
formly rotating solar-type stars. These models are com- 
plementary to the models we presented in Paper I for 
gas giant planets. In our previous analysis, we showed 
that the presence of a solid core in gas giants can lead to 
a rich spectrum of features in the frequency-dependence 
of Q' . Since there are extended ranges in the core size 
and spin periods of gas giant planets, the dispersion in 
their values of Q' may be large. Nevertheless, the typical 
magnitude of Q' for gas giants is comparable to that for 
solar-type stars. 

It is useful to make a parallel comparison between 
the theoretically determined and observationally inferred 
values of Q' . Owing to the technical challenges of the 
precision measurement of radial velocities, it is difficult 
to reduce the uncertainties in planetary orbital eccentric- 
ity to less than about 0.05. Many planets with periods 
less than 5 d have measured eccentricities less than 0.05. 
There are also supposedly eccentric planets within this 
period range. However, the apparent eccentricity mea- 
surement in such cases may be due to the perturbation 
of the host star's motion by additional planets in the 
system. 

Finally, we note that the brown dwarf companion 
HD 41004Bb has a measured eccentricity of 0.08 despite 
its period of only 1.33 d. The mass ratio of the system is 
about 0.05. The observed eccentricity can be preserved 
over the 1.6 Gyr age of the system provided both ob- 
jects have Q' > lO"^. The value Q'^ ^ 3 x 10^ required 
for the preservation of this brown dwarf implies that the 



M dwarf is unlikely to have been effective in damping 
the eccentricity of this system (see the discussions in 
Section [4.1.2p . The lower limit on Q'^ required for the 
brown dwarf not to damp the eccentricity is an order 
of magnitude above that derived for the Jovian planet 
considered in Paper I. The weak dissipation in this case 
may be due to the different internal structure of brown 
dwarfs. Hough waves cannot be excited in these fully 
convective bodies. Without the presence of a core, in- 
ertial waves take the form of smooth global modes and 
are excite d only in th e vicinity of particular resonant fre- 
quencies (lWull2005allbl l. 

5. SUMMARY AND DISCUSSION 

In this paper, we have calculated the excitation and 
dissipation of low-frequency tidal oscillations in uni- 
formly rotating solar-type stars and investigated the con- 
sequences for the circularization of binary stars and the 
orbital migration of close planetary companions. Our 
objectives were, first, to try to account for the efficiency 
of tidal dissipation inferred from the circularization pe- 
riods of binary stars and, second, to try to resolve the 
paradox that gas giant planets with very short periods 
can retained despite their intense tidal interaction with 
their host stars. 

The detailed analysis presented here is an extension of 
our previous investigation. In Paper I, we studied the 
excitation and dissipation of inertial waves and Hough 
waves in the convective interior and radiative surface 
layers of gas giant planets. In the present paper we 
have examined solar-type stars with convective envelopes 
and radiative interiors. Again, our calculation allows for 
two avenues of tidal dissipation. The disturbance excited 
in the convective zone, which takes the form of inertial 
waves if the tidal frequency is less than twice the spin fre- 
quency, is dissipated by turbulent viscosity. Hough waves 
are also excited at the interface between the radiative 
and convective regions and propagate towards the center 
of the star. We calculate the associated dissipation rate 
under the assumption that they do not reflect coherently 
from the center of the star. 

The radiative zone presents a reflecting barrier to in- 
ertial waves because the stable stratification strongly in- 
hibits radial motions of low frequency. In this way the 
radiative interior of the star has a similar role to the 
solid core in the planetary model. When the inertial 
waves propagate within a spherical annulus (rather than 
a full sphere) the tidal response is spatially intricate and 
enhanced dissipation occurs over extended ranges of fre- 
quency. The value of Q' can be reduced by up to four 
orders of magnitude when the Coriolis force is taken into 
account. 

The values of Q' derived for the binary circularization 
problem in the case of the present Sun with a spin period 
of 10 d (Fig. [3|) may still be an order of magnitude too 
large to account for the observed circularization periods 
of solar-type binary stars. However, the sensitivity of 
these results to the tidal and spin frequencies, the stellar 
model, and possible differential rotation means that the 
mechanisms we describe might be adequate to explain 
the observations when a detailed evolutionary study is 
carried out. 

The dependence of Q' on the tidal and spin frequencies 
potentially resolves an outstanding issue with regard to 
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the survival of gas giant planets with very short-period 
orbits around solar-type stars. These planets cannot sup- 
ply angular momentum through tidal transfer at a suf- 
ficient rate to compensate for the spin-down of the host 
stars due to magnetic braking. The spin frequencies of 
these stars are comparable to that of the Sun, which is at 
least an order of magnitude smaller than the tidal forc- 
ing frequency. Consequently inertial waves cannot be 
excited. In addition we estimate that the Hough waves 
excited by a hot Jupiter in the present Sun are unlikely 
to undergo nonlinear damping at the center of the star. 
If the Hough waves reflect coherently from the center and 
form global modes with narrow resonances, their contri- 
bution to orbital migration may be negligible. These very 
short-period planets can therefore be preserved despite 
their proximity to their host stars. Based on our models, 
we infer that gas giants with periods less than 1 d period 
may be preserved around solar-type stars. However, as 
conditions at the center of the star evolve, nonlinearity 
may set in at a critical age, resulting in a relatively rapid 
inspiral of the planet. 



It would be valuable to conduct numerical simulations 
of gravity (or Hough) waves approaching the center of 
a star to determine at what amplitude the reflection of 
the wave is affected by nonlinearity (cf. Appendix [X| . 
Unfortunately, currently available stellar models do not 
agree on the value of dN/dr at the center, which is of 
critical importance. 

Recently, a number of new planetary transit candidates 
were announced w ith periods less than 1 d in some cases 
(|Sahu et al.ll2006[ ). If confirmed, which seems difficult, 
these systems would present interesting constraints on 
theories of tidal evolution. 
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APPENDIX 

GRAVITY WAVES APPROACHING THE CENTER OF THE STAR 

Lo w-frequency gravity waves of small amplitude in a non-rotating star are described by the differential equation 
(e.g. iGoodman fc DicksonlHOQl 



d_ 
dr 



ld_ 
p dr 



{pr%) 



1) 1 



(Al) 



where — Re[^r(?')^/"(^, 0) exp(— iwt)] is the radial displacement, p(r) is the density and N(r) is the Brunt-Vaisala 
frequency. Very near the center of the star, the density and pressure are nearly uniform while g is linear in r. In this 
region p can be regarded as constant while N = Cr is linear in r. Let x = N/oj = {C/uj)r. Then the regular solution 
of the above equation for £ > 1 is 

OC x-^/^M^e{i+l)x), (A2) 

in terms of the Bessel function of order i' = ^+ 1- (The second solution, involving Y^, diverges as r ^ 0.) In particular, 
for £ = 2, 
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where A is a dimensionless amplitude. Since the unperturbed entropy profile is proportional to r^, the perturbed 
profile is proportional to — 2r^r- The wave overturns the stratification, according to linear theory, if the quantity 



r or 



(A4) 



The maximum value of this quantity occurs ai r — 0, 9 ~ 7r/2, and 2(j) — tot — 0, and is equal to A. Therefore A is 
a measure of the maximum nonlinearity of the wave. If A < 1 the entropy profile is never inver ted. [In this case the 
wave may still be disrupted by parametric instabilities, as noted bv IGoodman fc DicksonI lUll).] 
The energy flux in the incoming part of the wave is 

F=^^. (A5) 
72V6 C5 ^ ^ 

Equating this to the energy flux in gravity waves excited at the interface between the radiative and convective zones 
in a non-rotating solar model, we estimate 



A « 2500e 



P 



1/6 



(A6) 



for the largest eccentricity tide involved in the circularization of a binary with two synchronized 1 Mq stars, while 

for the tide excited in a s lowly rotating 1 Mp, s t ar by a planet of mass Mp and orbital period P. The first estimate 
is similar to that made bv IGoodman fc DicksonI ()1998f ). and implies that the waves always achieve nonlinearity in the 
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binary circularization problem if the eccentricity is measurable. (This remains true when the Coriolis force is taken 
into account.) The second estimate implies that the waves excited by hot Jupiters are of significant amplitude but 
probably do not achieve sufficient nonlinearity to disrupt the reflection. 
These results depend strongly on the stellar model. For stars of mass 0.5 Mq and age 5 Gyr we find 

/pX 1/6 

A^2.1el-\ (A8) 



for binary circularization and 



for hot Jupiters. Therefore even the binary circularization problem is unlikely to give rise to nonlinear tides if both 
components have 0.5 Mq. The main effect causing this difference is that A cx {dN/dr)c^^ and {dN/dr)c is about 
7 times larger in the present Sun than in the 0.5 Mq star. In general {dN/dr)c increases with stellar mass and age, 
but its value is not very accurately determined by current models of stellar evolution. 
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